Fractional vortices and composite domain walls in flat nanomagnets 
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We provide a simple explanation of complex magnetic patterns observed in ferromagnetic nanos- 
tructures. To this end we identify elementary topological defects in the field of magnetization: 
ordinary vortices in the bulk and vortices with half-integer winding numbers confined to the edge. 
Domain walls found in experiments and numerical simulations in strips and rings are composite 
objects containing two or more of the elementary defects. 



Topological defects 0, 0] greatly influence the proper- 
ties of materials by catalyzing or inhibiting the switch- 
ing between differently ordered states. In ferromagnetic 
nanoparticles of various shapes (e.g. strips Q and rings 
0), the switching process involves creation, propaga- 
tion, and annihilation of domain walls with complex in- 
ternal structure Here we show that these domain 
walls are composite objects made of two or more ele- 
mentary defects: vortices with integer winding numbers 
(n = ±1) and edge defects with fractional winding num- 
bers (n = ±1/2). The simplest domain walls are com- 
posed of two edge defects with opposite winding numbers. 
Creation and annihilation of the defects is constrained 
by conservation of a topological charge. This framework 
provides a basic understanding of the complex switching 
processes observed in ferromagnetic nanoparticles. 

In ferromagnets the competition between exchange and 
magnetic dipolar energies creates nonuniform patterns 
of magnetization in the ground state. Whereas the ex- 
change energy favors a state with uniform magnetization, 
magnetic dipolar interactions align the vector of magne- 
tization with the surface. In a large magnet a compro- 
mise is reached by the formation of uniformly magnetized 
domains separated by domain walls. In a nanomagnet 
magnetization varies continuously forming intricate yet 
highly reproducible textures, which include domain walls 
and vortices [E El 0] • Numerical simulations 0, |j| reveal 
a rich internal structure and complex dynamics of these 
objects. For example, collisions of two domain walls can 
have drastically different outcomes: complete annihila- 
tion or formation of other stable textures [j| . These puz- 
zling phenomena call for a theoretical explanation. 

An elementary picture of topological defects in nano- 
magnets with a planar geometry is presented in this Let- 
ter. It is suggested that the elementary defects are vor- 
tices with integer winding numbers and edge defects with 
half-integer winding numbers. All of the intricate tex- 
tures, including the domain walls, are composite objects 
made of two or more elementary defects. 

For simplicity, we consider a ferromagnet without 
intrinsic anisotropy, which is a good approximation 
for permalloy. The magnetic energy consists of two 
parts: the exchange contribution A J |Vm| d 3 r, where 
rh = M/|M| is the unit vector pointing in the direc- 



tion of magnetization M, and the magnetostatic energy 
(Ho/2) J | H | 2 d 3 r. The magnetic field H is related to the 
magnetization through Maxwell's equations, V x H = 
and V • (H + M) = 0. Apart from a few special cases 
(e.g. an ellipsoidal particle), finding configurations rii(r) 
of lowest energy is a difficult computational problem. 

Analytical treatment is nonetheless possible in a thin- 
film limit defined for a strip of width w and thick- 



ness t as 



t < w <C X 2 /t < wlog (w/t), 



(1) 

where A = yj A/fioM 2 is a magnetic length scale (5 nm 
in permalloy). Taking this limit yields three simplifica- 
tions: (a) magnetization lies in the plane of the film, 
m = (cos 9, sin 9, 0); (b) it depends on the in-plane coor- 
dinates x and y, but not on z; (c) the magnetic energy 
becomes a local functional of magnetization 0, [n| : 



E[m(r)]/At = [ |Vm| 2 
Jn 



d 2 r~(l/k) 



/ (m-r); 

JdQ 



dr. (2) 



The quantity E/At is a dimensionless energy, SI is the 
two-dimensional region of the film, dfl is its line bound- 
ary, r = (cos 9 T , sin 9 T , 0) is the unit vector parallel to 
the boundary, and A = 47rA 2 /i log (w/t) is the effective 
magnetic length in the thin-film geometry. The first — 
exchange — term in Eq. J5J is the familiar XY model 1] 
with the ground states obeying the Laplace equation 
V 2 = in the bulk. The second term expresses the 
magnetostatic energy due to magnetic charges at the film 
edge and sets the boundary conditions. 

Things are particularly simple in the limit A/w — > 0, 
when magnetization at the edge is forced to be paral- 
lel to the boundary, m = ±t. Solutions in this limit 
can be constructed by exploiting the analogy between 
the XY model and electrostatics in two dimensions Q, 
whereby the angle gradients are associated with com- 
ponents of the electric field, (E x ,E y ) = (d y 9,—d x 9). 
Topological defects — vortices with winding numbers n — 
±1 — become positive and negative point charges of unit 
strength. In the presence of an edge, solutions are ob- 
tained by the method of images. For instance, a single 
vortex (winding number n = ±1) located at (X, Y) in 
the infinite semiplane y > is 

6(x, y) = ± arctan ( — — ^ J ± arctan ( — — — ) . (3) 
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FIG. 1: (a) and (b) Edge defects with winding numbers n — +| and - 
— | edge defects, (d) A domain wall composed of two edge defects ±-| 
energy dominates. 



|. (c) A vortex (n = +1) offsets the influence of two 
(e) The — i defect in the limit where magnetostatic 



The image "charge" at (X, —Y) has the same sign as the 
original. Therefore, a vortex is repelled by the boundary. 

Another class of solutions [Fig. ^a) and (b)] has a 
singularity at the edge 



9(x, y) = ± arctan 



X 



(4) 



In the spirit of Eq. (|3J), this solution can be viewed as a 
halfvortex whose winding number ±5 is doubled by the 
reflection in the edge. Confinement to the edge occurs be- 
cause vortices with noninteger winding numbers cannot 
exist in the bulk. An attempt to peel a halfvortex away 
from the edge will make 9(x, y) a multivalued function: 



9(x, y) = ± — arctan 



Y 



X 



± — arctan 
2 



Y 



X 



(5) 

The multivaluedness means that the magnetization will 
be discontinuous across a line defect extending from the 
core of the halfvortex to the film boundary and thus pro- 
viding a potential V ex Y confining the halfvortex to the 
edge. 

The assignment of fractional winding numbers to edge 
defects can be justified in a more direct way: a superpo- 
sition of two edge defects and one vortex (n = +1) has 
zero total defect strength [Fig. Gfc)]. Hence n = — ^ for 
each of the edge defects. A model-independent, purely 
geometrical justification will be given below. 

In a strip \y\ < w/2, a domain wall interpolating be- 
tween the 9 = and ir ground states can be constructed 
out of two edge defects with opposite winding numbers 



±i and^ [Fig-Htd)]: 



tan#(x, y) = ± 



cos (iry/w) 



sinh (n(x — X)/w) 



(6) 



It is clear from the electrostatic analogy that the two 
defects experience an attractive "Coulomb" force (not 
strong enough to overcome the edge confinement). This 
force holds the composite domain wall together. Because 
the XY model itself contains no length scale, the extent 
of the domain wall along the strip is set by the width of 
the strip w. This is precisely the kind of a domain wall 
observed in numerical simulations 

The above solutions can be easily adopted to the cases 
when the ratio A/w is small but finite. The singular cores 
of the halfvortices reside outside the film, the distance A 
away from the edge flf)j |. 

The topological defects identified above survive well 
outside the range of applicability of the model J5J). Fig. [21 
shows simulated dynamics in a permalloy disk of the di- 
ameter w = 400 nm and thickness t = 20 nm. (We used 
the numerical package OOMMF Inequality Q fails 

badly, yet the bulk and edge defects are easily recognis- 
able and the winding numbers are clearly conserved. 

In wider and thicker strips (wt > A 2 ), the magneto- 
static energy breaks the symmetry between the defects 
with positive and negative winding numbers inherent to 
the XY model. In particular, the +5 edge defects have a 
higher magnetostatic energy than their — ^ counterparts. 
Evidence for this can be seen in numerical simulations 
showing a more substantial broadening of the cores in 
n 



•5 defects. 
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FIG. 2: Numerical simulation of magnetization dynamics in a permalloy disk. Color encodes the direction of the field rh (also 
shown by arrows). Convergence of different colors at a point signals the presence of a topological defect. In the first three 
panels the disk contains two +1 vortices and two — A edge defects. After the edge defects and a vortex mutually annihilate 
(last panel), a single vortex remains. 



In the magnetostatic regime wt ^> A 2 , the defects can 
again be constructed explicitly 01 . The magnetostatic 
energy is minimized if the density of magnetic charges 
vanishes in the bulk, V • m = 0, and on the surface, 
m = ±t. It can be checked that the +1 vortex in the 
bulk retains its shape [Fig. Cfc)] . The — \ edge defect 
looks like two domain walls emanating from a point at the 
edge in orthogonal directions [Fig.^e)]. The —1 vortex is 
a similar intersection of four domain walls. The V-shaped 
— h edge defects can be seen in numerical simulations . 
Domain walls in this limit are made of two — A defects 
and a +1 vortex between them [l3|. This type of the 
domain wall is favored when wt > CA 2 , where C « 130 is 
a rather large numerical parameter [|| weakly dependent 
on the thickness t. 

Although different in shape, the — \ edge defects in the 
exchange [Fig. ^b)] and magnetostatic [Fig. (He)] limits 



J 



have identical topological properties. Generally, as long 
as magnetization tends to align itself with the boundary, 
the edge defects are stable. They are manifested as kinks 
in magnetization rh: the latter rotates along the edge 
between the local tangential directions +r and — t. The 
winding number of a kink is defined as the line integral 
along the boundary: 



hi v < 

^ Jon 



6U • dr 



1 

±- 

2 



(7) 



for a single edge defect. Because the boundary dfl is a 
closed line, it always contains an even number of kinks. 
The sum of the winding numbers of edge defects is thus 
an integer related to the topological charge of bulk de- 
fects. In a simply connected (no holes) region O, 



edge If bulk 

Yn t = l-—d> V6-dr = l-— {d x d y 9 - d y d x B) d 2 r = 1 - V m 
t 27T J on m Jn i 



(8) 



r 



with integer rii for bulk defects. In a film with g holes, 

Y^m = i-g. (9) 

i 

The winding numbers are ±1 for bulk defects and ±A 
for edge defects. The half-integer value of the topolog- 
ical charge is directly related to the pres ence of a kink 
in magnetization at the edge 0, llil Il6l| . The halfvor- 
tices are analogs of the boojums thought to exist at an 
interface between the A and B phases of superfluid 3 He 

Bin. 

Conservation of topological charge © has important 
implications for the dynamics of magnetization in nano- 



magnets. For example, in rings of certain sizes the ground 
state contains no topological defects and has zero mag- 
netic dipole moment [Fig. Ola)] . The ring can be mag- 
netized by applying a strong in-plane magnetic field. 
Switching off the field leaves the ring in a metastable 
state with remnant magnetization containing two com- 
posite domain walls [Fig. [Sfb)] . By applying the mag- 
netic field in the opposite direction, the walls can be set 
in motion on a collision course and may annihilate, leav- 
ing the magnet in a ground state 0. However, direct 
annihilation of the two domain walls is impossible: both 
walls have the — \ defects at the inner edge of the ring 
and two such defects cannot annihilate on their own. Ac- 
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FIG. 3: Stable and metastable states of a magnetic nanor- 
ing in the thin-film limit, (a) One of the two ground states 
without topological defects, (b) State with a remnant magne- 
tization contains two composite domain walls with +i defects 
at the outer edge and — | defects at the inner edge. 



cordingly, in thin and narrow rings the two domain walls 
do not annihilate but instead form a "360° domain wall" 
7J. In thicker and wider rings, the annihilation does oc- 
cur and the ring returns to a ground state. To facilitate 
the annihilation, the edge defects in one of the domain 
walls must trade places. Because the edge defects them- 
selves cannot move into the bulk, they alter the signs by 
exchanging a vortex. The +^ defect emits a +1 vortex 
into the bulk and converts into a — \ defect jj^- The 
emitted vortex travels to the inner edge where it fuses 
with the — \ defect into a +i defect. Now the defects at 
both edges can annihilate directly. In thin and narrow 
rings this does not happen because vortex emission is for- 
bidden energetically |2j . Transient states with vortices in 
the bulk of the ring have been observed in simulations 
0. However, their role in catalysing the annihilation of 
domain walls has not been appreciated. 

We have uncovered a simple structure underlying intri- 
cate magnetic patterns in soft nanomagnets in a planar 
geometry. The patterns are formed by a few highly sta- 



ble topological defects: ordinary vortices in the bulk and 
halfvortices confined to the edge. The importance of edge 
defects in nanomagnets has been overlooked. Here we 
have demonstrated that the "transverse" domain walls 
found in thin and narrow strips are made of two such 
edge defects wiht opposite winding numbers. We will 
describe elsewhere ll3j a simple model of the "vortex" 
domain walls made of two — i edge defects and a +1 
vortex. It appears that the general approach taken in 
this paper may provide a basic model of complex magne- 
tization dynamics in nanomagnets by reducing it to the 
creation, propagation, and annihilation of a few topolog- 
ical defects. 
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